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Reorientation transition of ultrathin ferromagnetic films 
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We demonstrate that the reorientation transition from out-of-plane to in-plane magnetization 
with decreasing temperature as observed experimentally in Ni-films on Cu(OOl) can be explained 
on a microscopic basis. Using a combination of mean field theory and perturbation theory, we de- 
rive an analytic expression for the temperature dependent anisotropy. The reduced magnetization 
in the film surface at finite temperatures plays a crucial role for this transition as with increasing 
temperature the influence of the uniaxial anisotropies is reduced at the surface and is enhanced 
inside the film. 

PACS numbers: 68.35.Rh, 75.10.Hk, 75.30.Gw, 75.70.-i 



The direction of the magnetization of thin ferromag- 
netic films depends on various anisotropic energy contri- 
butions like,surface anisotropy fields which often favor an 
orientationEI perpendicular to the film, dipole interaction 
which favors an in-plane magnetization, and eventually 
anisotropy fields in the inner layers. As a consequence of 
these competing effects, a temperature driven reorienta- 
tion transition from an out-of-plane ordered state at low 
temperatures to an in-plane ordered state at high tem- 
peratures may be observed at appropriate chosen film 
thicknesses. Experimentally, this transition has 
studied in detail for various ultra-thin magnetic filmscl □. 
Recently, it was found by Schulz and BaberschkeB that 
ultra-thin Ni-films grown on Cu(OOl) show an opposite 
behavior: the magnetization is oriented in-plane for low 
temperatures and perpendicular at high temperatures. 

Phenomenological approaches for explaining the reori- 
entation transition usually start from the energy (or the 
free energy at finite temperatures) which is expanded 
in terms of the orientation of the magnetization vector 
relative to the film introducing temperature dependent 
anisotropy coefficients Ki(T). The temperature depen- 
dence of these coefficients is then studied experimentally 
(for a recent review seeQ). 

To better understand the mechanism responsible for 
the temperature driven transition, several investigations 
have been done in the framework of statistical spin mod- 
els. The advantage of this approach is that only a few 
microscopic parameters enter: besides an exchange inter- 
action the dipole interaction and an uniaxial anisotropy 
in the surface layers of the film. While Moschel et aZJll 
showed that the temperature dependence of the reorien- 
tation transition is well described qualitatively within a 
quantum mechanical mean field approach, most other au- 
thors focused on classical spin madels. Extended Monte 
Carlo simulations on mono-layeEa3a as well as pajean field 
calculations of both mono-layersEEl and bilayersO agree in 
the sense that a temperature driven reorientation transi- 
tion is obtained. Nevertheless, there is still a controvers, 
with respect to the order of this transition. While Chuil 
measured the expectation value of the components of the 
total magnetization and obtained a second order transi- 
tion for a monolayer we found, using an improved simula- 



tion algorithm and analyzing the Monte Carlo data with 
a histogram method, a transition of first order in agree- 
ment with the mean field calculations for this systemLl. 
Furthermore, we could show that the order of the transi- 
tion depends on the number of lavers and on the distri- 
bution of the uniaxial anisotropiestil. 

In all of these theoretical investigations a temperature 
driven reorientation transition from a out-of-plane state 
at low temperatures to an in-plane state at high temper- 
ature is found for appropriate sets of parameters, which 
is due to a competition of a positive surface anisotropy 
and the dipole interaction. TheL-interesting new result 
for ultra-thin Ni-films is arguedO to have its origin in 
a stress-induced uniaxial anisotropy energy in the in- 
ner layers with its easy axis perpendicular to the film. 
This anisotropy is in competition with the dipole in- 
teraction and a negative surface anisotropy. While the 
thickness-dependeat transition could be explained with 
these anisotropiea3, the origin of the more interesting 
temperature driven transition is not yet explained on a 
microscopic basis. Note that the_reversed reorientation 
recently found by Maclsaac et aZ.Ej has a different origin 
as it only occurs at vanishing exchange interaction. It 
is the purpose of this letter to show that the dipole in- 
teraction together with uniaxial anisotropies in the film 
indeed may lead to a temperature driven second order re- 
orientation transition from an in-plane magnetized film 
at low temperatures to a perpendicular magnetized film 
at high temperatures. 

After submission of this letter -pe became aware of a 
paper by Jensen and BennemannEj on the same topic. 
Starting from an expansion of the free energy in terms of 
uniaxial anisotropy and dipole interaction and employ- 
ing 1-h^n a mean field approximation following earlier 
worlo they calculated numerically the temperature of 
both types of reorientation transitions. In contrast to 
this calculation we develop in the present paper a fully 
selfconsistent mean field theory and analyse the reorien- 
tation transition within this approach since only within 
a nonlinear theory the canted phase and in particular its 
width can be analyzed. Additionally, a selfconsistent cal- 
culation of the quantities Ki{T) introduced below which 
are crucial for an understanding of the nature of the tran- 
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sition is not possible. A linearization of the free energy 
which is discussed in the last part of our paper agrees 
with the results of Ref.li3. Within such a linearized the- 
ory the approximate location of the transition can be 
obtained as a temperature somewhere within the canted 
phase but it is not possible to calculate the width of this 
phase. Thus, this calculation is only meaningful for a 
situation where the canted phase occupies a rather small 
temperature interval. Note, however, that this width can 
be quite large for instance for strongly varyina-anisotropy 
energies or certain parameter configurationalil in which 
case a nonlinear approach is necessary. 

The calculations are done in the framework of a classi- 
cal ferromagnetic Heisenberg model consisting of L two- 
dimensional layers on a simple cubic lattice. The Hamil- 
tonian reads 



n 



J 



- 3r„-''^(s*, • )(fy • Sj), (1) 



where Si = {sf , sf,sl) are spin vectors of unit length at 



position ri = [r 



in layer and fi 



n -fj. J 

is the nearest- neighbor exchange coupling constant, Dx 
is the uniaxial anisotropy which depends on the layer 
index A = 1 . . . L, and lo = ^qp? / Aira^ is the strength of 
the long range dipole interaction on a lattice with lat- 
tice constant a (/io is the magnetic permeability and /i is 
the effective magnetic moment of one spin) . All energies 
and temperatures are measured in units of J (fee — 1) 
which is fixed to J = 1 in this letter. Note that only sec- 
ond order uniaxial anisotropics D\ enter the Hamiltonian 
Eq. (U) . In our calculations we will restrict ourself to the 
case that all anisotropics are the same except at one sur- 
face, as this scenario is sufficient for explaining the basic 
physics of the temperature driven reorientation transi- 
tion. Furthermore we will focus on the case of L = 4 
layers. A systematic investigation of the parameter and 
thickness dependence of the reorientation transition and 
in particular a calcuLation of the corresponding phase di- 
agrams is under wayMI. 

In the following we assume translational invariance 
within the layers and therefore we set (s, ) — rri\ if Si is a 
spin in layer A. For the Hamiltonian Eq. (|l|) a molecular- 
field approximation is implemented resulting in L effec- 
tive one particle Hamiltonians from which the free energy 
functional can be obtained: the mean field in layer A is 
given by h\ = 'X.x^rfi^ where X^^ contains both ex- 
change and dipole interaction. With the order parameter 
M = {rfii, rriL) the Hamiltonian in layer A becomes 



njMF 



(M) 



z \2 



■mx - s A - Dx{s\) 



(2) 



Integrating this mean field Hamiltonian over the surface 
of the unit sphere in each layer yields the free energy per 
surface element. 




0.5 

Reduced Temperature 

FIG. 1. Magnetization components and anisotropy en- 
ergies for a Fe-type system with L = 4 layers. J = 1, 
Di/J = 14 X 10"^, Dx>i = 0, Lo/J = 38 X 10"^ The param- 
eters are based on Iron. 



^(T,M) = -r^log 

A=l 



ds 



■ (M)/T 



(3) 



Due to the in-plane rotational invariance of Eq. 

and thus the free energy in Eq. 



pends on the 2L components rrvi and m| of M and is 



can set vn\ 




stationary with respect to variations of these quantities. 
This variation is done in two steps: First we minimize the 
free energy Eq. with the constraint that the azimuth 
angle •& of the total magnetization m = 
fixed, and expand the resulting constrained free energy 
in powers of cos('!9) to give the angle-dependent free en- 
ergy 



-Kiir) cos*(??) - ... 



(4) 



with the reduced temperature r = T/Tc {Tc is the Curie 
temperature of the film) and temperature dependent ex- 
pansion coefficients Kiij). These quantities are usu- 
ally introduced phenomenologically. However in our ap- 
proach we can calculate these coefficients Ki (r) from the 
microscopic parameters of the system. The equilibrium 
free energy is then obtained as the minimum of Eq. (^) 
with respect to -d. 

In this notation, the two reorientation transition tem- 



peratures r^, where 
given by the conditions 



^2«) 



0, and T^, where ' 



0, are 

(5a) 
(5b) 



Fig. 1^ shows the temperature dependence of the compo- 
nents of the total magnetization m(r) and the anisotropy 
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FIG. 2. Magnetization components and anisotropy en- 
ergies for a Ni-type system with L = 4 layers. J = 1, 
DxjJ = -3.5 X 10-^ Dx>i/J = 1.5 x 10-^ cj/J = 5 x 10'^ 
The parameters are based on Nickel. 



Eq. ^ considering lj and Dx as small perturbations of 
the pure isotropic Heisenberg Hamiltonian which is jus- 
tified in view of the smallness of these parameters. We 
will only give the results of these calculations in this let- 
ter, the complete depvation will be reported in detail in 
a forthcoming paperEj. 

The total anisotropy K{t) of the system is defined as 
the difference of the free energies of the in-plane state 
and the out-of-plane state 



K{t)=T{t,7t/2)-T{t,0) 

^K2{T)+Ki{T) + ... 



(6) 



When we neglect the narrow canted phase, the reorienta- 
tion temperature is given by the condition if (r^) = 0. 
If the first derivative 9^_ftr(rr) < 0, we have a normal 
transition from out-of-plane to in-plane magnetization 
direction, otherwise the transition is reversed. 

In the framework of a perturbation theory we can 
derive an analytical expression for K{t) involving the 
absolute value of the layer magnetizations toa(t) and 
the fluctuations transversal to the magnetization direc- 
tion (7a(t) = ((sf both calculated with the unper- 
turbed Hamiltonian: 



coefBcients K2{t) and Ki{T) for a situation where one 
of the layers (the surface layer) has a positive uniaxial 
anisotropy, Di > 0, and the others are set to £'a>i = 0. 
Kq (t) and higher order terms are nearly two magnitudes 
smaller than K2{t) and therefore are not depicted. This 
is the situation encountered for instance in ultra-thin Fe- 
films. The ground-state magnetization of the system is 
perpendicular to the film. Increasing the temperature 
the magnetization switches continuously from the per- 
pendicular direction at a temperature to the in-plane 
direction at r^. 

For Ni-films on Cu(OOl) there is a positive uniaxial 
volume anisotropy in the inner layers favoring perpen- 
dicular orientation, and eventually a negative anisotropy 
on the surfaced. In competition to these energies is the 
dipole interaction which always favors in-plane magneti- 
zation. Fig. H shows the temperature dependence of the 
components of the total magnetization vector and the 
anisotropy coefficients Kiir) for a Ni-type system with 
L — A layers. The exchange interaction J is estimated 
from the Curie temperature of bulk Nickel, the dipole 
constant uj is calculated from the ground state magnetic 
moment and the lattice constant, and the anisotropy en- 
ergies are taken from the experiments. For these pa- 
rameters with increasing temperature the magnetization 
starts to cant at a temperature and reaches the per- 
pendicular state at as observed experimentally in Ni- 
films. These results were obtained numerically by solving 
the corresponding mean field equations. 

In order to understand both the normal and the re- 
versed reorientation transition we additionally applied 
a perturbation theory to the mean field Hamiltonian 



K{t) ^ Kq{T) + K^{t) 

L 
A=l 

3a; ^ 

~T ^ mx{T)^\x-y\my{T). (7) 

A,A' = 1 

The constants $a contain the effective dipole inter- 
action between the layers and can be calculated nu- 
merically to give $0 = 9.0336, $i = -0.32750, and 
ci>,>i = 0(e-2-(^-i)). 

At the critical temperature K{t) vanishes and hence 
K[t) must be curved in order to have another zero at a 
temperature < 1. Furthermore, a positive curvature 
is necessary for a normal reorientation transition while a 
negative curvature of K{t) is necessary for a reversed re- 
orientation. Hence we will focus on the second derivatives 
of Eq. (|^) and start with the dipole part Km (t) . It turns 
out that d^Kmij) > for all film thicknesses and tem- 
peratures since uj is positive and the main contribution of 
the sum is proportional to X]a"^a(''") which always has 
a negative curvature. Thus the dipole interaction always 
favors the normal reorientation and can never lead to a 
reversed transition in an exchange dominated system. 

Now we will examine Kq{T). First note that q\{Q) — 
and q\{l) — 1/3 in the unperturbated case. For L — 1 
and L = 2 we have q\{T) — r/i in mean field approxi- 
mation, and then Kq{T) — J2\ D\{1 — t), i.e. the second 
derivative vanishes in this case. Consequently in systems 
with L < 2 layers we only find normal reorientation tran- 
sitions from out-of-plane direction at low temperatures to 
an in-plane direction at higher temperatures. 
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this case only a reorientation transition from out-of-plane 
to in-plane can occur, provided the exchange interaction 
is large with respect to the uniaxial anisotropics and the 
dipole interaction. When the film thickness L > 2, the 
magnetization is not homogeneous through the film as 
the surface layers have a reduced magnetization at fi- 
nite temperatures. This leads to an enhancement of the 
transversal fluctuations at the surface and to a reduc- 
tion of these fluctuations in the inner of the film. Hence 
the influence of the uniaxial anisotropics is reduced at 
the surface and enhanced inside the film favoring a spin 
orientation parallel to the easy axis of the inner layer. 
This effect may lead to a temperature driven reorienta- 
tion transition of the type observed in Nickel. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft through Sondcrforschungsbereich 166. The 
authors would like to thank S. Liibeck and A. Moschel 
for fruitful discussions. 



FIG. 3. Total anisotropy K{t) and its two parts Kq[T) 
and Km(T) from Eq. ^ for a Ni-type system. The model 
parameters are the same as in Fig. S. 



This is not the case for L > 2 layers since then, due to 
the reduced surface magnetization at finite temperatures, 
the transversal fluctuations in the surface layers qs (r) are 
enhanced (qs(r) > r/3, d^Qsir) < 0). Combined with a 
negative surface anisotropy this may lead to a neg- 
ative curvature of Kqij). Furthermore, the transversal 
fluctuations in the inner layers are reduced by this effect 
i<lb{T) < t/3, d^qbir) > 0) and, when combined with a 
positive uniaxial anisotropy in the inner layers enhance 
the negative curvature of Kg (r) . 

In Fig. H K{t) is depicted together with the two com- 
peting parts Kqij) and Kmij) from Eq. (Q) for the same 
parameters as in Fig. ^. A transition is obtained with in- 
creasing temperature because Kqlr) tends slower to zero 
than Kmir). 

In summary we have shown that the temperature 
driven reorientation transitions seen in ultra-thin ferro- 
magnetic films are well described within a mean field 
approximation if second order uniaxial anisotropics and 
the dipole interaction are included in the Hamiltonian. 
In particular we can relate the unusual transition seen 
in Ni-films to a microscopic model in which a positive 
uniaxial anisotropy energy is present in the inner lay- 
ers. Additionally we can calculate the parameters Ki(T) 
usually introduced phenomenologically from microscopic 
parameters of the system. 

The L = 4 layer film considered serves as a simple 
system showing the Ni-type transition, while the Fe-type 
transition is already observed in mono layers if the pa- 
rameters are adjusted properly. This has a rather inter- 
esting physical origin: For systems with L — 1 or L — 2 
layers the unperturbed system is homogeneous as every 
lattice site has the same environment. It turns out that in 
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